## [1] "explicated variable of regression : canopy_cover"
## [1] "for Guinean_forest-savanna_regression_canopy_cover.RDS"
## [2] "for Northern_Congolian_Forest-Savanna_regression_canopy_cover.RDS"
## [3] "for Sahelian_Acacia_savanna_regression_canopy_cover.RDS"
## [4] "for Southern_Congolian_forest-savanna_regression_canopy_cover.RDS"
## [5] "for West_Sudanian_savanna_regression_canopy_cover.RDS"
## [6] "for Western_Congolian_forest-savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## [1] "########################################"
## [1] "########################################"
## [1] "below, stancode for Guinean_forest-savanna_regression_canopy_cover.RDS"
## // generated with brms 2.20.4
## functions {
## /* zero-inflated beta log-PDF of a single response
## * Args:
## * y: the response value
## * mu: mean parameter of the beta distribution
## * phi: precision parameter of the beta distribution
## * zi: zero-inflation probability
## * Returns:
## * a scalar to be added to the log posterior
## */
## real zero_inflated_beta_lpdf(real y, real mu, real phi, real zi) {
## row_vector[2] shape = [mu * phi, (1 - mu) * phi];
## if (y == 0) {
## return bernoulli_lpmf(1 | zi);
## } else {
## return bernoulli_lpmf(0 | zi) +
## beta_lpdf(y | shape[1], shape[2]);
## }
## }
## /* zero-inflated beta log-PDF of a single response
## * logit parameterization of the zero-inflation part
## * Args:
## * y: the response value
## * mu: mean parameter of the beta distribution
## * phi: precision parameter of the beta distribution
## * zi: linear predictor for zero-inflation part
## * Returns:
## * a scalar to be added to the log posterior
## */
## real zero_inflated_beta_logit_lpdf(real y, real mu, real phi, real zi) {
## row_vector[2] shape = [mu * phi, (1 - mu) * phi];
## if (y == 0) {
## return bernoulli_logit_lpmf(1 | zi);
## } else {
## return bernoulli_logit_lpmf(0 | zi) +
## beta_lpdf(y | shape[1], shape[2]);
## }
## }
## // zero-inflated beta log-CCDF and log-CDF functions
## real zero_inflated_beta_lccdf(real y, real mu, real phi, real zi) {
## row_vector[2] shape = [mu * phi, (1 - mu) * phi];
## return bernoulli_lpmf(0 | zi) + beta_lccdf(y | shape[1], shape[2]);
## }
## real zero_inflated_beta_lcdf(real y, real mu, real phi, real zi) {
## return log1m_exp(zero_inflated_beta_lccdf(y | mu, phi, zi));
## }
## }
## data {
## int<lower=1> N; // total number of observations
## vector[N] Y; // response variable
## int<lower=1> K; // number of population-level effects
## matrix[N, K] X; // population-level design matrix
## int<lower=1> Kc; // number of population-level effects after centering
## int prior_only; // should the likelihood be ignored?
## }
## transformed data {
## matrix[N, Kc] Xc; // centered version of X without an intercept
## vector[Kc] means_X; // column means of X before centering
## for (i in 2:K) {
## means_X[i - 1] = mean(X[, i]);
## Xc[, i - 1] = X[, i] - means_X[i - 1];
## }
## }
## parameters {
## vector[Kc] b; // regression coefficients
## real Intercept; // temporary intercept for centered predictors
## real<lower=0> phi; // precision parameter
## real<lower=0,upper=1> zi; // zero-inflation probability
## }
## transformed parameters {
## real lprior = 0; // prior contributions to the log posterior
## lprior += student_t_lpdf(Intercept | 3, 0, 2.5);
## lprior += gamma_lpdf(phi | 0.01, 0.01);
## lprior += beta_lpdf(zi | 1, 1);
## }
## model {
## // likelihood including constants
## if (!prior_only) {
## // initialize linear predictor term
## vector[N] mu = rep_vector(0.0, N);
## mu += Intercept + Xc * b;
## mu = inv_logit(mu);
## for (n in 1:N) {
## target += zero_inflated_beta_lpdf(Y[n] | mu[n], phi, zi);
## }
## }
## // priors including constants
## target += lprior;
## }
## generated quantities {
## // actual population-level intercept
## real b_Intercept = Intercept - dot_product(means_X, b);
## }
## [1] "########################################"
## [1] "########################################"
## [1] "########################################"
## [1] " "
## [1] " "
## [1] "Guinean_forest-savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## Family: zero_inflated_beta
## Links: mu = logit; phi = identity; zi = identity
## Formula: canopy_cover ~ fire_freq_std + mean_precip_std
## Data: table_region (Number of observations: 1725)
## Draws: 3 chains, each with iter = 10000; warmup = 2000; thin = 10;
## total post-warmup draws = 2400
##
## Population-Level Effects:
## Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## Intercept -1.95 0.06 -2.06 -1.84 1.00 2329 2393
## fire_freq_std -0.03 0.02 -0.07 0.02 1.00 2512 2410
## mean_precip_std 0.12 0.03 0.07 0.17 1.00 2246 2185
##
## Family Specific Parameters:
## Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## phi 5.28 0.19 4.93 5.68 1.00 2300 2519
## zi 0.11 0.01 0.10 0.13 1.00 2349 2292
##
## Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
## and Tail_ESS are effective sample size measures, and Rhat is the potential
## scale reduction factor on split chains (at convergence, Rhat = 1).
## [1] "########################################"
## [1] "Northern_Congolian_Forest-Savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## Family: zero_inflated_beta
## Links: mu = logit; phi = identity; zi = identity
## Formula: canopy_cover ~ fire_freq_std + mean_precip_std
## Data: table_region (Number of observations: 243)
## Draws: 3 chains, each with iter = 10000; warmup = 2000; thin = 10;
## total post-warmup draws = 2400
##
## Population-Level Effects:
## Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## Intercept -1.77 0.36 -2.49 -1.06 1.00 2062 2093
## fire_freq_std -0.00 0.06 -0.13 0.11 1.00 2320 2300
## mean_precip_std 0.10 0.17 -0.24 0.41 1.00 2021 2090
##
## Family Specific Parameters:
## Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## phi 6.72 0.62 5.57 7.98 1.00 2172 2174
## zi 0.06 0.02 0.03 0.09 1.00 2462 2175
##
## Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
## and Tail_ESS are effective sample size measures, and Rhat is the potential
## scale reduction factor on split chains (at convergence, Rhat = 1).
## [1] "########################################"
## [1] "Sahelian_Acacia_savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## Family: zero_inflated_beta
## Links: mu = logit; phi = identity; zi = identity
## Formula: canopy_cover ~ fire_freq_std + mean_precip_std
## Data: table_region (Number of observations: 5563)
## Draws: 3 chains, each with iter = 10000; warmup = 2000; thin = 10;
## total post-warmup draws = 2400
##
## Population-Level Effects:
## Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## Intercept -3.89 0.04 -3.96 -3.81 1.00 2460 2389
## fire_freq_std 0.18 0.03 0.12 0.23 1.00 2266 2215
## mean_precip_std 0.85 0.06 0.73 0.97 1.00 2353 2288
##
## Family Specific Parameters:
## Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## phi 27.54 1.24 25.14 29.96 1.00 2436 2332
## zi 0.68 0.01 0.67 0.69 1.00 2484 2090
##
## Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
## and Tail_ESS are effective sample size measures, and Rhat is the potential
## scale reduction factor on split chains (at convergence, Rhat = 1).
## [1] "########################################"
## [1] "Southern_Congolian_forest-savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## Family: zero_inflated_beta
## Links: mu = logit; phi = identity; zi = identity
## Formula: canopy_cover ~ fire_freq_std + mean_precip_std
## Data: table_region (Number of observations: 47)
## Draws: 3 chains, each with iter = 10000; warmup = 2000; thin = 10;
## total post-warmup draws = 2400
##
## Population-Level Effects:
## Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## Intercept -1.68 2.43 -6.50 3.17 1.00 2565 2352
## fire_freq_std -0.04 0.15 -0.34 0.25 1.00 2875 2352
## mean_precip_std -0.20 1.18 -2.55 2.05 1.00 2539 2409
##
## Family Specific Parameters:
## Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## phi 4.47 1.51 2.01 7.78 1.00 2263 2328
## zi 0.51 0.07 0.38 0.64 1.00 2388 2116
##
## Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
## and Tail_ESS are effective sample size measures, and Rhat is the potential
## scale reduction factor on split chains (at convergence, Rhat = 1).
## [1] "########################################"
## [1] "West_Sudanian_savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## Family: zero_inflated_beta
## Links: mu = logit; phi = identity; zi = identity
## Formula: canopy_cover ~ fire_freq_std + mean_precip_std
## Data: table_region (Number of observations: 3277)
## Draws: 3 chains, each with iter = 10000; warmup = 2000; thin = 10;
## total post-warmup draws = 2400
##
## Population-Level Effects:
## Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## Intercept -3.52 0.05 -3.61 -3.42 1.00 2378 2285
## fire_freq_std 0.08 0.01 0.05 0.11 1.00 2402 2369
## mean_precip_std 0.90 0.04 0.82 0.98 1.00 2507 2256
##
## Family Specific Parameters:
## Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## phi 9.97 0.32 9.36 10.63 1.00 2315 2251
## zi 0.23 0.01 0.21 0.24 1.00 2245 2288
##
## Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
## and Tail_ESS are effective sample size measures, and Rhat is the potential
## scale reduction factor on split chains (at convergence, Rhat = 1).
## [1] "########################################"
## [1] "Western_Congolian_forest-savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## Family: zero_inflated_beta
## Links: mu = logit; phi = identity; zi = identity
## Formula: canopy_cover ~ fire_freq_std + mean_precip_std
## Data: table_region (Number of observations: 259)
## Draws: 3 chains, each with iter = 10000; warmup = 2000; thin = 10;
## total post-warmup draws = 2400
##
## Population-Level Effects:
## Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## Intercept -1.91 0.33 -2.54 -1.25 1.00 2474 2412
## fire_freq_std -0.13 0.05 -0.23 -0.03 1.00 2465 2251
## mean_precip_std 0.20 0.21 -0.22 0.58 1.00 2423 2348
##
## Family Specific Parameters:
## Estimate Est.Error l-95% CI u-95% CI Rhat Bulk_ESS Tail_ESS
## phi 3.51 0.39 2.79 4.30 1.00 2551 2455
## zi 0.22 0.03 0.17 0.28 1.00 2147 2145
##
## Draws were sampled using sampling(NUTS). For each parameter, Bulk_ESS
## and Tail_ESS are effective sample size measures, and Rhat is the potential
## scale reduction factor on split chains (at convergence, Rhat = 1).
## [1] "########################################"
## [1] "Gamma regressions for Guinean_forest-savanna_regression_rh98.RDS"
## [2] "Gamma regressions for Northern_Congolian_Forest-Savanna_regression_rh98.RDS"
## [3] "Gamma regressions for Sahelian_Acacia_savanna_regression_rh98.RDS"
## [4] "Gamma regressions for Southern_Congolian_forest-savanna_regression_rh98.RDS"
## [5] "Gamma regressions for West_Sudanian_savanna_regression_rh98.RDS"
## [6] "Gamma regressions for Western_Congolian_forest-savanna_regression_rh98.RDS"
## [1] "########################################"
## [1] "########################################"
## [1] "########################################"
## [1] "Guinean_forest-savanna_regression_rh98.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1] 1725 3200
## [1] "(nb_donnes I * nb_iter_mcmc J)"
## [1] "mean(table_region$rh98)"
## [1] 8.17064
## [1] "sd(table_region$rh98)"
## [1] 4.228092
## [1] "mean(simulations[,j]) ( truth = 8.171 )"
## [1] 8.061
## [1] "sd(simulations[,j]) ( truth = 4.228 )"
## [1] 4.058
## [1] "mean(simulations[,j]) ( truth = 8.171 )"
## [1] 8.1
## [1] "sd(simulations[,j]) ( truth = 4.228 )"
## [1] 4.143
## [1] "mean(simulations[,j]) ( truth = 8.171 )"
## [1] 8.231
## [1] "sd(simulations[,j]) ( truth = 4.228 )"
## [1] 4.174
## [1] "mean(simulations[,j]) ( truth = 8.171 )"
## [1] 8.561
## [1] "sd(simulations[,j]) ( truth = 4.228 )"
## [1] 4.492
## [1] "mean(simulations[,j]) ( truth = 8.171 )"
## [1] 8.023
## [1] "sd(simulations[,j]) ( truth = 4.228 )"
## [1] 4.201
## [1] "Northern_Congolian_Forest-Savanna_regression_rh98.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1] 243 3200
## [1] "(nb_donnes I * nb_iter_mcmc J)"
## [1] "mean(table_region$rh98)"
## [1] 8.251706
## [1] "sd(table_region$rh98)"
## [1] 3.52155
## [1] "mean(simulations[,j]) ( truth = 8.252 )"
## [1] 8.045
## [1] "sd(simulations[,j]) ( truth = 3.522 )"
## [1] 3.504
## [1] "mean(simulations[,j]) ( truth = 8.252 )"
## [1] 7.762
## [1] "sd(simulations[,j]) ( truth = 3.522 )"
## [1] 3.474
## [1] "mean(simulations[,j]) ( truth = 8.252 )"
## [1] 8.521
## [1] "sd(simulations[,j]) ( truth = 3.522 )"
## [1] 3.923
## [1] "mean(simulations[,j]) ( truth = 8.252 )"
## [1] 7.937
## [1] "sd(simulations[,j]) ( truth = 3.522 )"
## [1] 3.17
## [1] "mean(simulations[,j]) ( truth = 8.252 )"
## [1] 8.253
## [1] "sd(simulations[,j]) ( truth = 3.522 )"
## [1] 3.75
## [1] "Sahelian_Acacia_savanna_regression_rh98.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1] 5563 3200
## [1] "(nb_donnes I * nb_iter_mcmc J)"
## [1] "mean(table_region$rh98)"
## [1] 3.049537
## [1] "sd(table_region$rh98)"
## [1] 1.187586
## [1] "mean(simulations[,j]) ( truth = 3.05 )"
## [1] 3.057
## [1] "sd(simulations[,j]) ( truth = 1.188 )"
## [1] 0.963
## [1] "mean(simulations[,j]) ( truth = 3.05 )"
## [1] 3.037
## [1] "sd(simulations[,j]) ( truth = 1.188 )"
## [1] 0.972
## [1] "mean(simulations[,j]) ( truth = 3.05 )"
## [1] 3.056
## [1] "sd(simulations[,j]) ( truth = 1.188 )"
## [1] 1.017
## [1] "mean(simulations[,j]) ( truth = 3.05 )"
## [1] 3.057
## [1] "sd(simulations[,j]) ( truth = 1.188 )"
## [1] 0.999
## [1] "mean(simulations[,j]) ( truth = 3.05 )"
## [1] 3.076
## [1] "sd(simulations[,j]) ( truth = 1.188 )"
## [1] 1.035
## [1] "Southern_Congolian_forest-savanna_regression_rh98.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1] 47 3200
## [1] "(nb_donnes I * nb_iter_mcmc J)"
## [1] "mean(table_region$rh98)"
## [1] 3.938866
## [1] "sd(table_region$rh98)"
## [1] 2.086246
## [1] "mean(simulations[,j]) ( truth = 3.939 )"
## [1] 4.539
## [1] "sd(simulations[,j]) ( truth = 2.086 )"
## [1] 2.147
## [1] "mean(simulations[,j]) ( truth = 3.939 )"
## [1] 4.376
## [1] "sd(simulations[,j]) ( truth = 2.086 )"
## [1] 2.142
## [1] "mean(simulations[,j]) ( truth = 3.939 )"
## [1] 3.702
## [1] "sd(simulations[,j]) ( truth = 2.086 )"
## [1] 0.973
## [1] "mean(simulations[,j]) ( truth = 3.939 )"
## [1] 3.446
## [1] "sd(simulations[,j]) ( truth = 2.086 )"
## [1] 1.511
## [1] "mean(simulations[,j]) ( truth = 3.939 )"
## [1] 5.172
## [1] "sd(simulations[,j]) ( truth = 2.086 )"
## [1] 3.096
## [1] "West_Sudanian_savanna_regression_rh98.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1] 3277 3200
## [1] "(nb_donnes I * nb_iter_mcmc J)"
## [1] "mean(table_region$rh98)"
## [1] 5.728088
## [1] "sd(table_region$rh98)"
## [1] 3.15554
## [1] "mean(simulations[,j]) ( truth = 5.728 )"
## [1] 5.838
## [1] "sd(simulations[,j]) ( truth = 3.156 )"
## [1] 3.069
## [1] "mean(simulations[,j]) ( truth = 5.728 )"
## [1] 5.853
## [1] "sd(simulations[,j]) ( truth = 3.156 )"
## [1] 3.239
## [1] "mean(simulations[,j]) ( truth = 5.728 )"
## [1] 5.647
## [1] "sd(simulations[,j]) ( truth = 3.156 )"
## [1] 2.967
## [1] "mean(simulations[,j]) ( truth = 5.728 )"
## [1] 5.742
## [1] "sd(simulations[,j]) ( truth = 3.156 )"
## [1] 3.096
## [1] "mean(simulations[,j]) ( truth = 5.728 )"
## [1] 5.674
## [1] "sd(simulations[,j]) ( truth = 3.156 )"
## [1] 3.021
## [1] "Western_Congolian_forest-savanna_regression_rh98.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1] 259 3200
## [1] "(nb_donnes I * nb_iter_mcmc J)"
## [1] "mean(table_region$rh98)"
## [1] 6.137749
## [1] "sd(table_region$rh98)"
## [1] 4.607916
## [1] "mean(simulations[,j]) ( truth = 6.138 )"
## [1] 6.332
## [1] "sd(simulations[,j]) ( truth = 4.608 )"
## [1] 3.828
## [1] "mean(simulations[,j]) ( truth = 6.138 )"
## [1] 6.284
## [1] "sd(simulations[,j]) ( truth = 4.608 )"
## [1] 4.183
## [1] "mean(simulations[,j]) ( truth = 6.138 )"
## [1] 6.204
## [1] "sd(simulations[,j]) ( truth = 4.608 )"
## [1] 3.881
## [1] "mean(simulations[,j]) ( truth = 6.138 )"
## [1] 6.442
## [1] "sd(simulations[,j]) ( truth = 4.608 )"
## [1] 3.943
## [1] "mean(simulations[,j]) ( truth = 6.138 )"
## [1] 5.531
## [1] "sd(simulations[,j]) ( truth = 4.608 )"
## [1] 3.084
## [1] "Beta regressions for Guinean_forest-savanna_regression_canopy_cover.RDS"
## [2] "Beta regressions for Northern_Congolian_Forest-Savanna_regression_canopy_cover.RDS"
## [3] "Beta regressions for Sahelian_Acacia_savanna_regression_canopy_cover.RDS"
## [4] "Beta regressions for Southern_Congolian_forest-savanna_regression_canopy_cover.RDS"
## [5] "Beta regressions for West_Sudanian_savanna_regression_canopy_cover.RDS"
## [6] "Beta regressions for Western_Congolian_forest-savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## [1] "########################################"
## [1] "########################################"
## [1] "Guinean_forest-savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1] 1725 2400
## [1] "(nb_donnes I * nb_iter_mcmc J)"
## Le chargement a nécessité le package : gtools
##
## Attachement du package : 'gtools'
## Les objets suivants sont masqués depuis 'package:brms':
##
## ddirichlet, rdirichlet
## [1] "mean(table_region$canopy_cover)"
## [1] 0.1314814
## [1] "sd(table_region$canopy_cover)"
## [1] 0.14466
## [1] "mean(simulations[,j]) ( truth = 0.131 )"
## [1] 0.162
## [1] "sd(simulations[,j]) ( truth = 0.145 )"
## [1] 0.15
## [1] "mean(simulations[,j]) ( truth = 0.131 )"
## [1] 0.168
## [1] "sd(simulations[,j]) ( truth = 0.145 )"
## [1] 0.152
## [1] "mean(simulations[,j]) ( truth = 0.131 )"
## [1] 0.167
## [1] "sd(simulations[,j]) ( truth = 0.145 )"
## [1] 0.158
## [1] "mean(simulations[,j]) ( truth = 0.131 )"
## [1] 0.159
## [1] "sd(simulations[,j]) ( truth = 0.145 )"
## [1] 0.148
## [1] "mean(simulations[,j]) ( truth = 0.131 )"
## [1] 0.153
## [1] "sd(simulations[,j]) ( truth = 0.145 )"
## [1] 0.144
## [1] "Northern_Congolian_Forest-Savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1] 243 2400
## [1] "(nb_donnes I * nb_iter_mcmc J)"
## [1] "mean(table_region$canopy_cover)"
## [1] 0.1649144
## [1] "sd(table_region$canopy_cover)"
## [1] 0.1334132
## [1] "mean(simulations[,j]) ( truth = 0.165 )"
## [1] 0.166
## [1] "sd(simulations[,j]) ( truth = 0.133 )"
## [1] 0.143
## [1] "mean(simulations[,j]) ( truth = 0.165 )"
## [1] 0.199
## [1] "sd(simulations[,j]) ( truth = 0.133 )"
## [1] 0.158
## [1] "mean(simulations[,j]) ( truth = 0.165 )"
## [1] 0.169
## [1] "sd(simulations[,j]) ( truth = 0.133 )"
## [1] 0.138
## [1] "mean(simulations[,j]) ( truth = 0.165 )"
## [1] 0.158
## [1] "sd(simulations[,j]) ( truth = 0.133 )"
## [1] 0.125
## [1] "mean(simulations[,j]) ( truth = 0.165 )"
## [1] 0.161
## [1] "sd(simulations[,j]) ( truth = 0.133 )"
## [1] 0.127
## [1] "Sahelian_Acacia_savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1] 5563 2400
## [1] "(nb_donnes I * nb_iter_mcmc J)"
## [1] "mean(table_region$canopy_cover)"
## [1] 0.005161202
## [1] "sd(table_region$canopy_cover)"
## [1] 0.02202149
## [1] "mean(simulations[,j]) ( truth = 0.005 )"
## [1] 0.464
## [1] "sd(simulations[,j]) ( truth = 0.022 )"
## [1] 0.319
## [1] "mean(simulations[,j]) ( truth = 0.005 )"
## [1] 0.478
## [1] "sd(simulations[,j]) ( truth = 0.022 )"
## [1] 0.317
## [1] "mean(simulations[,j]) ( truth = 0.005 )"
## [1] 0.477
## [1] "sd(simulations[,j]) ( truth = 0.022 )"
## [1] 0.314
## [1] "mean(simulations[,j]) ( truth = 0.005 )"
## [1] 0.472
## [1] "sd(simulations[,j]) ( truth = 0.022 )"
## [1] 0.317
## [1] "mean(simulations[,j]) ( truth = 0.005 )"
## [1] 0.485
## [1] "sd(simulations[,j]) ( truth = 0.022 )"
## [1] 0.315
## [1] "Southern_Congolian_forest-savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1] 47 2400
## [1] "(nb_donnes I * nb_iter_mcmc J)"
## [1] "mean(table_region$canopy_cover)"
## [1] 0.04984602
## [1] "sd(table_region$canopy_cover)"
## [1] 0.1063901
## [1] "mean(simulations[,j]) ( truth = 0.05 )"
## [1] 0.393
## [1] "sd(simulations[,j]) ( truth = 0.106 )"
## [1] 0.256
## [1] "mean(simulations[,j]) ( truth = 0.05 )"
## [1] 0.339
## [1] "sd(simulations[,j]) ( truth = 0.106 )"
## [1] 0.29
## [1] "mean(simulations[,j]) ( truth = 0.05 )"
## [1] 0.435
## [1] "sd(simulations[,j]) ( truth = 0.106 )"
## [1] 0.283
## [1] "mean(simulations[,j]) ( truth = 0.05 )"
## [1] 0.378
## [1] "sd(simulations[,j]) ( truth = 0.106 )"
## [1] 0.305
## [1] "mean(simulations[,j]) ( truth = 0.05 )"
## [1] 0.358
## [1] "sd(simulations[,j]) ( truth = 0.106 )"
## [1] 0.297
## [1] "West_Sudanian_savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1] 3277 2400
## [1] "(nb_donnes I * nb_iter_mcmc J)"
## [1] "mean(table_region$canopy_cover)"
## [1] 0.05300694
## [1] "sd(table_region$canopy_cover)"
## [1] 0.08641946
## [1] "mean(simulations[,j]) ( truth = 0.053 )"
## [1] 0.119
## [1] "sd(simulations[,j]) ( truth = 0.086 )"
## [1] 0.13
## [1] "mean(simulations[,j]) ( truth = 0.053 )"
## [1] 0.117
## [1] "sd(simulations[,j]) ( truth = 0.086 )"
## [1] 0.126
## [1] "mean(simulations[,j]) ( truth = 0.053 )"
## [1] 0.11
## [1] "sd(simulations[,j]) ( truth = 0.086 )"
## [1] 0.123
## [1] "mean(simulations[,j]) ( truth = 0.053 )"
## [1] 0.115
## [1] "sd(simulations[,j]) ( truth = 0.086 )"
## [1] 0.127
## [1] "mean(simulations[,j]) ( truth = 0.053 )"
## [1] 0.117
## [1] "sd(simulations[,j]) ( truth = 0.086 )"
## [1] 0.129
## [1] "Western_Congolian_forest-savanna_regression_canopy_cover.RDS"
## [1] "########################################"
## [1] "dim(linear_predictors_for_one_beta_draw_per_column)"
## [1] 259 2400
## [1] "(nb_donnes I * nb_iter_mcmc J)"
## [1] "mean(table_region$canopy_cover)"
## [1] 0.1079882
## [1] "sd(table_region$canopy_cover)"
## [1] 0.1649516
## [1] "mean(simulations[,j]) ( truth = 0.108 )"
## [1] 0.168
## [1] "sd(simulations[,j]) ( truth = 0.165 )"
## [1] 0.19
## [1] "mean(simulations[,j]) ( truth = 0.108 )"
## [1] 0.19
## [1] "sd(simulations[,j]) ( truth = 0.165 )"
## [1] 0.182
## [1] "mean(simulations[,j]) ( truth = 0.108 )"
## [1] 0.187
## [1] "sd(simulations[,j]) ( truth = 0.165 )"
## [1] 0.185
## [1] "mean(simulations[,j]) ( truth = 0.108 )"
## [1] 0.179
## [1] "sd(simulations[,j]) ( truth = 0.165 )"
## [1] 0.195
## [1] "mean(simulations[,j]) ( truth = 0.108 )"
## [1] 0.187
## [1] "sd(simulations[,j]) ( truth = 0.165 )"
## [1] 0.192
